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In conventional treatments, predictions from fixed-order perturbative QCD calculations cannot 
be fixed with certainty due to ambiguities in the choice of the renormalization scale as well as the 
renormalization scheme. In this paper we present a general discussion of the constraints of the 
renormalization group (RG) invariance on the choice of the renormalization scale. We adopt the 
RG based equations, which incorporate the scheme parameters, for a general exposition of RG 
invariance, since they simultaneously express the invariance of physical observables under both 
the variation of the renormalization scale and the renormalization scheme parameters. We then 
discuss the self-consistency requirements of the RG, such as reflexivity, symmetry, and transitivity, 
which must be satisfied by the scale-setting method. The Principle of Minimal Sensitivity (PMS) 
requires the slope of the approximant of an observable to vanish at the renormalization point. This 
criterion provides a scheme-independent estimation, but it violates the symmetry and transitivity 
properties of the RG and does not reproduce the Gell-Mann-Low scale for QED observables. The 
Principle of Maximum Conformality (PMC) satisfies all of the deductions of the RG invariance - 
reflectivity, symmetry, and transitivity. Using the PMC, all non-conformal {/3^}-terms (1Z stands 
for an arbitrary renormalization scheme) in the perturbative expansion series are summed into 
the running coupling, and one obtains a unique, scale-fixed, scheme-independent prediction at any 
finite order. The PMC scales and the resulting finite-order PMC predictions are both to high 
accuracy independent of the choice of initial renormalization scale, consistent with RG invariance. 
Moreover, after PMC scale-setting, the residual initial scale-dependence at fixed order due to 
unknown higher-order {/3i}-terms can be substantially suppressed. The PMC thus eliminates a 
serious systematic scale error in pQCD predictions, greatly improving the precision of tests of the 
Standard Model and the sensitivity to new physics at collider and other experiments. 
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I. INTRODUCTION 

Given the perturbative series for a physical quantity 

n 

p n = C a?(M)+^C i (M)< +i (M) (1) 

expanded to n-th order in the QCD strong coupling con- 
stant a s (/i), the renormalization scale /x must be speci- 
fied in order to obtain a definite prediction. The com- 
mon practice adopted in the literature is to simply guess 
a renormalization scale fj, = Q, keep it fixed during the 
calculation (Q is usually assumed to be a typical momen- 
tum transfer of the process), and then vary it over an 
arbitrary range, e.g. [Q/2, 2Q], in order to ascertain the 
scale-uncertainty. However there are many weak points 
of this conventional scale-setting method : 

1. Although the infinite perturbative series p n ^oo 



summed to all orders is renormalization-scale in- 
dependent, the scale dependence from a s (jj,) and 
Ci(fj) do not exactly cancel at finite order, leading 
to a renormalization scale ambiguity. 

2. The fixed-order estimate in the conventional proce- 
dure is also scheme dependent; i.e., different choice 
of renormalization scheme 1Z will lead to different 
theoretical estimates. This is the well-known renor- 
malization scheme ambiguity [lMTTj. 

3. The conventional scale choice can give unphysical 
results: For example, for the case of W^-boson plus 
three-jet production at the hadronic colliders, tak- 
ing fj, to be the T^-boson transverse energy, the con- 
ventional scale-setting method even predicts nega- 
tive QCD cross-section at the next-to-leading-order 
(NLO) [Hill. 

4. As has been shown in Ref. [l4j |. taking an incorrect 
renormalization scale underestimates the top quark 
forward-backward asymmetry at the Tevatron. 
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5. It should be recalled that there is no ambiguity in 
setting the renormalization scale in QED. 
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In QED, the coupling a(q 2 ) is convention- 
ally denned in the Gell Mann-Low (GM- 
L) scheme [l5[ from the potential between 
heavy charges, and it is normalized at q 2 = 
to the fine-structure constant a(0) ~ 
1/137.0359... [Tf3]. 

Due to the Ward-Takahashi identity [17], the 
divergences in the vertex and fermion wave- 
function corrections cancel, and the ultraviolet 
divergence associated with the vacuum polar- 
ization defines a natural scale for the coupling 
constant a{q 2 ) with q 2 being the squared mo- 
mentum transfer for the photon propagator. 
This fact was first observed by Gell-Mann and 
Low [l5j : i.e., in the standard GM-L scheme, 
the renormalization scale is simply the virtu- 
ally of the exchanged photon. 
For example, the renormalization scale for the 
electron-muon elastic scattering based on one- 
photon exchange is the virtuality of the ex- 



changed photon, i.e. ^gm- 



t 



One 



can of course choose any initial renormaliza- 
tion scale to for calculating the QED ampli- 
tude; however, the final result will not depend 
on the choice of to, since 



a{t) 



a(t ) 



where 



n(t,t ) 



[i-n(t,t ) 



[nffl - u(t Q )} 
[i - n(to)] 



(2) 



naturally sums all vacuum polarization con- 
tributions, both proper and improper, to the 
dressed photon propagator. (Here LI(t) = 
II(i, 0) is the sum of proper vacuum polariza- 
tion insertions, subtracted at t = 0.) The 
invariance of the result on the initial scale 
to is the property used to derive the Callan- 
Symanzik equations [HI, [l9j]. There is, there- 
fore, no reason to vary /iqm-l by a factor of 
1/2 or 2, since the photon virtuality t is the 
unique, optimized scale in the GM-L scheme. 

• The renormalization scale in QED is unique 
in any scheme including dimensional regular- 
ization; different schemes can be connected to 
the GM-L scheme by commensurate scale rela- 
tions (CSRs) , a topic which we discuss be- 
low. The computation of higher-order 
functions is thus important for perturbative 
calculations at the highest orders [2H - I23] ]. 

6. There are uncancelled large logarithms, as well 
as "renormalon" terms in higher orders which di- 
verge as [n!(/3^) n a™] [H[. The convergence of 
the perturbative scries is thus problematic using 
conventional scale-setting. For certain processes 



such as the top-quark pair production, it is found 
that the total cross-section for the (gg)-channel, 
qq — > t + t, at the next-to-next-to-leading order 
(NNLO) is about 50% of the NLO cross-section us- 
ing conventional scale-setting [25T427T ] . Thus, to de- 
rive a dependable perturbative estimate, one evi- 
dently needs to do even higher order calculations. 

7. The conventional estimate shows a strong depen- 
dence on the choice of the renormalization scale 
[i. It is clearly artificial to guess a renormalization 
scale /1 = Q and to study its uncertainty by simply 
varying /1 in the arbitrary range [Q/2,2Q]. Why 
is the scale uncertainty estimated only by varying 
a factor of 1/2 or 2, and not, say, 10 times Q ? For 
example, Ref. 28] shows that after including the 
first and second order corrections to several deep 
inelastic sum rules which are due to heavy flavor 
contributions, it is found that the effective scale 
[i ~ 6.5 Q, where the typical scale Q — m with m 
being the corresponding heavy quark mass. 

Moreover, sometimes, there are several choices for 
the typical momentum transfer of the process, all 
of which can be taken as the renormalization scale, 
such as the heavy-quark mass, the collision energy 
of the subprocess, etc. Which invariant provides 
the correct theoretical estimate ? 

Using conventional scale-setting, there is no def- 
inite answer to these questions. One may argue 
that the correct renormalization scale for the fixed- 
order prediction can be decided by comparing with 
the experimental data, but this surely is process de- 
pendent and greatly depresses the predictive power 
of the pQCD theory. 

Thus, in summary, the conventional scale-setting as- 
signs an arbitrary range and an arbitrary systematic error 
to fixed-order pQCD predictions. In fact, as we discuss 
in this article, this ad hoc assignment of the range and 
associated systematic error is unnecessary and, in fact, 
can be eliminated. 

One may ask: For a general fixed-order calculation, 
what is the correct "physical" scale or optimized scale ? 
To our understanding, it should provide a prediction in- 
dependent of the renormalization scheme and the choice 
of initial scale. A pictorial representation of what is 
the optimized renormalization scale is shown in Fig.tjTJ, 
where the electron-muon elastic scattering through one- 
photon exchange is taken as an illustration. In the GM-L 
scheme, the optimized scale [j, 2 pt = t which corresponds 
to the scale- invariant value a(t). Moreover, by using the 
proper scale-setting method, such as the newl y su ggested 
Principle of Maximum conformality (PMC) [25J, |29I431| , 
the prediction is also scheme independent and the argu- 
ment of the coupling in different schemes have the correct 
displacement. For example, by using the PMC procedure 
for QED one obtains the correct displacement between 
the argument of the coupling in the MS scheme relative 
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FIG. 1: Pictorial representation of the optimized renormal- 
ization scale jLt op t- Taking electron-muon elastic scattering 
through one-photon exchange as an example: In the GM-L 
scheme, the optimized scale is fj^pt = t which corresponds to 
the scale-invariant value a(t). As a comparison, the values of 
a at fixed-orders for different choice of to = Qi (i = 1, 2, 3, ■ • •) 
are shown by thin-and-solid curves. 



to the GM-L scheme at one loop 



ocgm-lK 



(t) = a^e-^H) . 



(3) 



As a comparison, the values of a at fixed-orders for 
different choice of to = Qf (i = 1, 2, 3, • • ■) are shown by 
thin-and-solid curves in Fig.(JTJ. A particular choice of to 
using conventional scale-setting may lead to a value of a 
close to a(t), but this would only be a lucky guess and not 
the correct answer. As one includes higher-and-higher 
orders, the guessed scale will lead to a better estimate. 
In fact, when doing the perturbative calculation up to 
infinite order, any choice of to will lead to the correct 
value a(t) as required by the RG invariance. However, 
if one chooses to = t, the complete all-orders result is 
obtained from the onset. 

Does there exist such an optimized renormalization 
scale for a general process in non-Abelian QCD ? If it 
does exist, how can one set it at finite order in a system- 
atic and process-independent way ? This is not an easy 
task. Various scale-setting procedures have been pro- 
posed since the 1980's for deriving an optimized scale, 
such as Fastest Apparent Convergence (FAC) Ji™ 2 , 
the Principle of Minimum Sensitivity (PMS) 0-[l3> tne 
Brodsky-Lepage-Mackenzie (BLM) [llj procedure and its 
extended versions such as the dressed skeleton expan- 
sion J34|, 35 1, the sequential se-BLM and x-BLM meth- 
ods [MlJ etc., and the PMC. A short review of FAC, 
BLM and PMS can be found in Ref. 39]. In principle, 
the correctness of a scale-setting method can be judged 
by the experimental data. However, as we shall dis- 
cuss, there are self-consistency theoretical requirements 



1 The displacement for higher-order corrections can be obtained 
by carefully dealing with the differences of the {/3i}-series under 
different renormalization schemes. 

2 As argued by Grunberg [f| and Krasnikov 33] , it is better to be 
called as the RG-improved effective coupling method. For sim- 
plicity, we retain the name as FAC as suggested by Stevenson Q . 



which shed light on the reliability of the scale-setting 
method 

Clearly, the prediction for any physical observable 
must be independent of the choice of renormalization 
scheme; this is the central pro p ert y of the renormaliza- 
tion group (RG) invariance [la, fl E). [41H 43|. As we shall 
discuss, the RG based equations [7H10| which incorporate 
the scheme parameters provide a convenient way for es- 
timating both the scale- and scheme- dependence of the 
QCD predictions for a physical process 0413, US Ell- 
in this paper, we will utilize such RG based equations 
for a general discussion of the RG-invariance. We will 
discuss in detail the self-consistency requirements of the 
RG [13], such as refiexivity, symmetry and transitivity, 
which must be satisfied by a scale-setting method. We 
will then show whether the scale-setting methods, FAC, 
BLM/PMC and PMS, satisfy these requirements. 

The remaining parts of this paper are organized as fol- 
lows: in Sec. II, we give a general demonstration of the 
RG-invariance with the help of the RG based equations. 
In Sec. Ill, we discuss the self-consistency requirements 
for a scale-setting method, where a graphical explana- 
tion of these requirements is also given. In Sec. IV and 
Sec. V, we present a detailed discussion on PMC and 
PMS scale-setting methods, respectively. Sec. VI pro- 
vides a summary. 



II. RENORMALIZATION-GROUP BASED 
EQUATIONS AND THE RENORMALIZATION 
GROUP INVARIANCE 

The scale dependence of the running coupling in gauge 
theory is controlled by the RG equation 







d In [J? 



47T 



i=0 



00 



4tt 



i+2 



(4) 



where the superscript 1Z stands for an arbitrary renormal- 
ization scheme, such as MS scheme [45[ , MS scheme [13] , 
MOM schem e ||[ , etc.. Note that the /^-functions for 
the MS and MS schemes are the same [47j | . Various 
terms in 0o~> P'i'i ' ' correspond to one- loop and two- 
loop • • • contributions respectively. In general, the 
are scheme-dependent and depend on the quark mass rrif. 
According to the decoupling theorem, a quark with mass 
rrif fi can be ignored, and we can often neglect m/- 
terms when rrif <C t 1 - Then, for every renormalization 
scale /x, one can divide the quarks into active ones with 
rrij = and inactive ones that can be ignored. Within 
this framework, it is well-known that the first two co- 
efficients /3q\ are universal; i.e., (3^ = H ~ 2n//3 and 
0^ = 102 — 38n//3 for n/-active flavors. Hereafter, we 
simply write them as fto and /3i . It is noted that an ana- 
lytic extension of ctjf s which incorporates the finite-mass 
quark threshold effects into the running of the coupling 
has been suggested in Ref. (48|. However, numerically, it 
is found that taking finite quark mass effects into account 
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analytically in the running, rather than using a fixed nt 
between thresholds, leads to effects of the order of one 
percent for the one- loop running coupling [4^. Here we 
will work with the conventional {/3^}-functions. 

It will be convenient to use the first two universal co- 
efficients Po and f3i to rescale the coupling constant and 
the scale-parameter in Eq.Q. By rescaling the coupling 
constant and the scale parameters as [44( 



Pi 



af and t-r 



4tt/3 



Pi 
Pi 



In// , 



one can express the RG equation Q into a simpler canon- 
ical form 



da n 



-{a n f [l + a n + cf(a n ) 2 + cf(a n f + 



(5) 

where cf = PfP^IP\ for i = 2, 3, ■ • •. 

As an extension of the ordinary coupling constant, one 
can define a universal coupling constant a(r, {ci}) to in- 
clude the dependence on the scheme parameters {ct}, 
which satisfies the following extended RG based equa- 
tions [3 



da 

P(a, {ct}) = — = -a 2 [l + a + c 2 a 2 
or 



c 3 a 



da f a x n+ dx 

Pn(a, {ci}) = — = -P(a, {a}) / 



(6) 



(7) 



The scale-equation ^ can be used to evolve the universal 
coupling function from one scale to another. The scheme- 
equation (|7J| , which was first suggested by Stevenson Q , 
can be used to relate the coupling functions under differ- 
ent schemes by changing {ci}. A solution of the scale- 
equation up to four- loop level has been given in Ref. [30| , 
which agrees with that of the conventional RG-equation 
obtained in the literature, cf. Ref. [49(. By comparing 
Eq.© with Eq.©, there exists a value of r = tr for 
which 



(r n ) = a(m,{cf}). 



This shows that any coupling constant a n (vji) can be 
expressed by the universal coupling constant a(r, {c,}) 
under proper correspondence; i.e. the coupling constant 
a n (T-ji) can be treated as a special case of the univer- 
sal coupling constant a(r, {q}): Any usual coupling con- 
stant a n (T-n) is equal to a universal coupling a(riz, {ci}) 
by setting {d} to be {cf-}, since both coupling constants 
satisfy the same RG equation by using the same scheme 
parameters. 

Grunberg has pointed out that 4 6] any pcrturbativcly 
calculable physical quantity can be used to define an ef- 
fective coupling constant by incorporating the entire ra- 
diative corrections into its definition. The effective cou- 
pling constant satisfies the same RG equation (and hence 
the same RG based equations) as the usual (universal) 
coupling constant. Thus, the running behavior for both 
the effective coupling constant and the usual (universal) 
coupling constant are the same if their RG based equa- 
tions are calculated under the same choice of scheme pa- 
rameters. This idea has later been discussed in detail by 
Refs. [10, [HI- Such an effective coupling constant can be 
used as a reference to define the renormalization proce- 
dure, such as MS scheme, MS scheme, etc.. 

The RG-invariance states that a physical quantity 
should be independent of the renormalization scale and 
renormalization scheme [IH [ltl I4H - I43I ] . This shows that 
if the effective coupling constant ci(tr, {cf-}) corresponds 
to a physical observable, it should be independent of any 
other scale t$ and any other scheme parameters {cf}; 
i.e. 



8t s 

da(TR,{cf}) 
def 



[scale invariance] , (9) 
[scheme invariance] . (10) 



(8) 



Demonstration: We provide an intuitive demonstration 
for the RG invariance from the above RG based equa- 
tions. Given two effective coupling constants (/(t-r., {cf}) 
and a(Tg,{cf}) defined under two different schemes TZ 
and S, one can expand a(r-ji, {cf}) in a power series of 
a(rs, {cf}) through a Taylor expansion: 



a{Tn,{ c i }) = a(r s +T, {c t + a}) 



+ 



2! 



cPa\ 
<9t 2 ) 



?2 + 2 



/ 8 2 a 
\ drdci 



TCi 



•E 



2 a 
deidej J s 1 



1 

3! 



d 3 a 



,(11) 



where r — t-r — t$, Ci = cf — cf and the subscript S next to the partial derivatives means they are evaluated 
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at the point (rs, {cf }). 

The right-hand side of Eq. fTTj) can be regrouped ac- 
cording to the different orders of scheme-parameters {c{\. 



where n stands for the highest perturbative order for a 
fixed-order calculation. It is noted that Eq. ([T2|) can be 
further simplified with the help of RG equations (|6l7p . 
If we set ti — y oo, the right-hand-side of Eg. (TT2"]) tends 
to zero, and we obtain the scale- invariance equation (|9|). 
This shows that if a{rn, {cf-}) corresponds to a physi- 
cal observable (corresponding to the case of n — > oo), 
it will be independent of any other scale T5. Similarly, 
doing the first derivative of aiju, {cf-}) with respect to 
the scheme-parameter cj , one can obtain the scheme- 
invariance equation (fTUj) . 

In another words, if one uses an effective coupling con- 
stant a(rs, {cf}) under the renormalization scheme S 
and with an initial renormalization scale {t$} to pre- 
dict the value of a(Tn, {cf}), the RG-invariances (|9ll0p 
tell us that 

• if we have summed all types of cf -terms (or equiv- 
alcntly the {/?f}-terms) into the coupling con- 
stant, as is the case of an infinite-order calcula- 
tion, then our final prediction of a{r-R,,{cf}) will 
be independent of any choice of initial scale ts and 



a{m,{cf}) = 



Using the RG scale-equation (|6]), the right- hand- 
side of the above equation can be rewritten as per- 
turbative series of a(rg, {cf}), whose coefficient at 
each order is a {/3^}-series. 



If one considers N c to be an analytic variable, then the 
scale-setting known from the non-Abelian theory SU(N C ) 
must agree with the Abelian QED theory at N c — > 0. 
This shows that above discussions are also suitable for 
QED; i.e. by taking the limit iV c — > at fixed a = Cpa s 
with Cf = (N„ — 1)/2N C , we effectively return to the 
QED case [53. l53j. 



After differentiating both side of Eq. (fTTj) over Tg, we ob- 
tain 



(12) 

I 

renormalization-scheme S. 

• According to Eq. (fT2|) . for a fixed-order estimation 
(i.e. n ^ 00), there is some residual initial-scale de- 
pendence. This is reasonable: as shown by Eg . (fTTj) . 
for a fixed-order calculation, the unknown {f3f}- 
terms in the higher orders are necessary to cancel 
the scale dependence from the lower-order terms. 

If we can find a proper way to sum up all the 
known- type of {/3f }-terms into the coupling con- 
stant, and at the same time suppressing the contri- 
butions from those unknown- type of {/3f }-terms ef- 
fectively, such residual initial scale dependence can 
be greatly suppressed. The PMC has been designed 
for such purpose [29l - [3lj ]. whose properties will be 
discussed in more detail in the following sections. 

• If setting all the differences of the renormalization 
scheme parameters, c\ = (i — 1,2, •••), Eq. fTTj) 
returns to a scale-expansion series for the coupling 
constant expanding over itself but specified at an- 
other scale; i.e. 



I 

III. SELF-CONSISTENCY REQUIREMENTS 
FOR A SCALE-SETTING METHOD 



It has been noted that if one knows how to set the 
optimal scale in all cases, then one can translate the re- 
sult freely from one scheme to another scheme through 
proper scale relations I34L 1 3 51 . This observation has later 
been emphasized in Ref. [20| , where the scale transforma- 
tion among different schemes are called "commensurate 
scale relations" (CSRs). It shows that even though the 
expansion coefficients could be different under different 
renormalization schemes, after a proper scale-setting, one 
can find a relation between the effective renormalization 
scales which ensures the total result remain the same 



da(Tn,{cf}) 9("+ 1 )a(T S ,{cf})f" 



dr s 



(n+l) 



^ d^a(r s ,{cf}) f 



'dr 



(n) 



(n-1)! 



I 

(w , {cf } ) + («2>) , + ' (^p) ,° + 1 p*p) ^ + ■ ■ ■ • as) 
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under any renormalization schemes. For simplicity, fol- 
lowing the suggestion of Ref . , we omit the scheme pa- 
rameters in the coupling constant in discussing the self- 
consistent requirements for a scale-setting method, but 
will retrieve them when necessary. 

In principle, the correctness of a scale-setting method 
can be judged by experimental data. However, it has 
been suggested that some self-consistency requirements 
can shed light on the reliability of the scale-setting 
method [40| . in which some initial discussions have been 
presented. These requirements together with their expla- 
nations are listed in the following: 

1. Existence and Uniqueness of the renormaliza- 
tion scale [i. Any scale-setting method must satisfy 
these two requirements. 

2. Reflexivity. Given an effective coupling a s (/i) 
specified at a renormalization scale /i, we can ex- 
press it in terms of itself but specified at another 
renormalization scale fi' , 



a a {n) = a s {n') + fi')al(fi') + 



(14) 



where /i') oc ln(/t 2 //x 2 ). When the scale // 
is chosen to be /x, the above equation reduces to a 
trivial identity. 

From the scale- invariance ©, up to infinite orders, 
we have 



(9 In// 2 



= 0. 



(15) 



This, inversely, means that if a s (/x) is known (say, 
a experimentally measured effective coupling), and 
we try to use the above perturbative equation to 
"predict" a s (p) from itself, then any deviation of 
// from n would lead to an inaccurate result due to 
the truncation of expansion series. 

More explicitly, for a fixed-order expansion with 
the highest perturbative-order n, from Ea. (TT2l . we 
obtain 



da. 



,(/i) (W/m' 2 ) 



<91n/Lt' 2 



<9(bV 2 )(™ 



+i) ' 



of 



This shows, generally, the right-hand-side 
Eq. ([T4|) depends on \i' at any fixed-order. 

Thus, to get a correct fixed-order estimate for 
a s (/i), a self-consistency scale-setting must take 
the unique value p! = fj, on the right-hand-side of 
Eq. (fT4)) . If a scale-setting satisfies such property, 
we say it is reflexive. 

It is found that the Reflexivity is a basic require- 
ment for a self-consistency scale-setting method 



and for the physical (effective) coupling constant 
a s (n), which provides the necessary condition 
for the following two properties Symmetry and 
Transitivity; i.e. if a scale-setting does not satisfy 
the Reflexivity, it cannot satisfy the following two 
properties Symmetry and Transitivity either. 

3. Symmetry. Given two different effective coupling 
constants a s i(/ii) and a s2 (/z 2 ) under two different 
renormalization schemes, we can expand any one of 
them in terms of the other: 



0^2(^2) 



a S 2(m) + r 12(^1,^2)04(^2) + 
a s i(/ii) + r 2 i(^2,Mi)asi(Mi) + 



After a general scale-setting, we have 



a s i(pi) 
0^2(^2) 

Note that, 



"32(^2) +^12(^1,^2)^2(^2/ 



• The new effective scales \i\ 2 may or may not 
be equal to fXi t 2, depending on the choice of 
the scale-setting method. The coefficients f i 2 
and ?i\ are changed accordingly in order to 
obtain a consistent result. 

• We have implicitly set the effective scales at 
NLO-level to be equal to the LO ones. We will 
adopt this choice throughout the paper. The 
effective scales for the highest-order terms are 
usually taken as the same effective scales at 
the one-lower-order, since they are the scales 
strictly set by the known-terms [30L HH . 

Setting = A 2 i^i and p\ = \\2P2, if 

Ai 2 A 21 = 1 , (16) 

we say that the scale-setting is symmetric. 

Explanation: 

If /i 2 = A 2 i/ii and [i\ = Ai 2 /i 2 , we obtain 

= a s2 (X2i^i) + fi 2 (ni, A 2 i^i)a 2 2 (A 2 i^i) H (17) 

and 
"52(^2) 

= a s i(Ai 2 ^2) + hiifJ-2, Ai 2 ^ 2 )a 2 !(Ai 2 ^ 2 ) H (18) 

As a combination of Eqs. (|17ll8|) . we obtain 
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a a i(pi) = a s i(Ai 2 A 2 iMi) + [hsiMi, A21 Mi) + r 2 i(A 2 i^i, Ai 2 A 2 iMi)] asi(^i2A2iMi) + 

I 



(19) 



From the Reflexivity property, if a scale-setting is 
symmetric, i.e. satisfying Eq. dTB)) . we will obtain 



ri2(^l,A* 2 ) + ^2l(yU2,Mi) = 0, 



(20) 



and vice versa. This shows that the Symmetry 
property (Til))) and the relation (|2U|) are mutually 
necessary and sufficient conditions. 

The Symmetry feature is necessary since it fur- 
ther gives us a unique relation for the scales before 
and after the scale-setting, 

= M1M2 • 

4. Transitivity. Given three effective coupling con- 
stants a s i(fii), a s2 (/i 2 ), and a s3 (/x 3 ) under three 
renormalization schemes, we can expand any one 
of them in terms of the other; i.e. 

a s i(ni) = a s2 (^ 2 ) +ri 2 (/ii,/i2)as 2 (M2) H , 

a s2 (/i 2 ) = a s3 (fi 3 ) + r 23 (^ 2 ,^ 3 )o4(/j 3 ) H , 

"33(^3) = a s i(Mi) + r 3 i(/i 3 ,/ii)asi(/xi) H . 

After a general scale-setting, we obtain 
a,i(Mi) = a S 2(/i 2 ) + ri 2 (^i,^ 2 )a^ 2 (^ 2 ) H , 



0^2(^2) = a S 3{n* 3 ) + f 23(^2,^3)^3(^3) H > 

a S 3(M3) = a a i(Mi) + ^13(^3,^1)^1(^1) H • 

Setting = A 2 i/ii, M3 = -^32^2 and /u| = Ai 3 ^ 3 , if 

Ai3A3 2 A 2 i = 1 . (21) 

we say that the scale-setting is transitive. 

Explanation: 

If ^ 2 = A 2 i/Ui, ^3 = A 32 ^ 2 and fi\ = Ai 3 ^ 3 , we 
obtain 

= Q. s2 (\2iHi) +ri 2 (/ii,A 2 i/xi)Q!g 2 (A2i/xi) H ,(22) 

0^2(^2) 

= a S 3(A 3 2M2) + r 23 (^ 2 , A3 2 /i 2 )a^3(A 32 M 2 ) H ,(23) 

0^3(^3) 

= a sl (X 13 fi 3 ) + f 31 0*3, Ai 3 /i 3 )a^ 1 (Ai 3 /i 3 ) H .(24) 

As a combination of Eqs. (l23l23l24[) . we obtain 



«si(mi) = a s i(Ai 3 A 32 A 2 i/ii) + a^ 1 (Ai 3 A 32 A 2 iMi) x 

[^3i(A 32 A 2 i/ii, Ai 3 A 32 A 2 i^i) + f 23 (A 2 i/xi, A 32 A 2 ipii) + fi 2 (/xi, A 2 i/ii)] + 

I 



(25) 



From the Reflexivity property, if a scale-setting is 
transitive, i.e. satisfying Eq. (|2"Tj) . we will obtain 



h2{Hl,l4) + f 2 3 0*2) M3) + ? 3l(M3i/ i l) = °, ( 26 ) 

and vice versa. This shows that the Transitivity 
property ([2~Tj) and the relation (|2^|) are mutually 
necessary and sufficient conditions. 

The Transitivity property shows that under a 
proper scale-setting method, we have A21 = A 23 A 3 i, 
which means that the scale ratio A 2 i for any two 
effective couplings a s \ and a s2 is independent of 
the choice of an intermediate effective coupling a s3 
under any renormalization scheme. Thus the re- 
lation between any two observables is independent 
of the choice of renormalization scheme. In fact, 
the Transitivity property provides the theoretical 



foundation for the existence of CSRs among differ- 
ent physical observables [20l |. 

The Transitivity feature gives us a unique relation 
for all the scales before and after the scale-setting, 



M1M2M3 = M1M2M3 ■ 

The Transitivity property is very important for 
a self-consistency scale-setting, which is a natural 
requirement from the RG invariance. It has al- 
ready been pointed out that why the renormaliza- 
tion group is called a "group" is mainly because of 
such Transitivity property (4lTj43j . 

The Transitivity property (l2~Tj) can be extended 
to an arbitrary number of effective coupling con- 
stants; i.e. if we have nth effective coupling con- 
stants which are related with similar manner as 
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above, then their transitivity relation is 

AinA„(„_i) • • • A32A21 = 1. (27) 

One may observe that the Symmetry is a spe- 
cial case of Transitivity, since if setting a S 3(//3) = 
0^1(^1); we have An = 1 and fn(^X)A*i) = due 
to the Renexivity, which thus changes the tran- 
sitive relation A13A32A21 = 1 into the symmetric 
relation A12A21 = 1. 

As a summary, a scale-setting method that satis- 
fies Existence and Uniqueness of the renormalization 
scale, Renexivity, Symmetry, and Transitivity effec- 
tively establishes equivalent relations among all the effec- 
tive coupling constants, and thus, among all the physical 
observables. 



A. A graphic explanation of these requirements 

In this subsection, we present a more intuitive expla- 
nation of these requirements based on the universal cou- 
pling a(r, {cj}) and the RG based Eqs. (|6l7[) . 

In the RG based equations (|6I7I) . there is no explicit 
reference to the QCD parameters, such as the number of 
colors or the number of active-flavors. Therefore, aside 
from its infinite dimensional character, a(r, {c^}) is just a 
mathematical function like, say, Bessel functions or any 
other special functions [44|. In practice, due to the un- 
known higher order scheme parameters {q}, wc need to 
truncate the beta function /3(a, {q}) and solve the uni- 
versal coupling constant a(r, {c^}) in a finite-dimensional 
subspace; i.e. we need to evaluate a(r, {c^}) in a subspace 
where higher order {ci}-terms are zero. In principle, this 
function can be computed to arbitrary degree of preci- 
sion, limited only by the truncation of the fundamental 
beta function. 

In this formalism, any two effective coupling constants 
can be related by some evolution path on the hyper- 
surface defined by a(r,{ci}). In Fig.© we illustrate 
the paths which represent the operations of Renexiv- 
ity, Symmetry and Transitivity. We can pictorially 
visualize that the evolution paths satisfy all these three 
self-consistency properties. A closed path starting and 
ending at the same point A represents the operation of 
identity. Since the predicted value does not depend on 
the chosen path, if the effective coupling constant at A 
is a a , after completing the path we will also end up with 
an effective coupling a a- Similarly, if we evolve as at 
5 to a value ac at C, we are guaranteed that when we 
evolve ac at C back to the point B, the result will be as- 
Hence, the evolution equations also satisfy Symmetry. 
Transitivity follows in a similar manner; i.e. going di- 
rectly from D to E gives the same result as going from 
D to E through a third point F. 

In the following two sections, we will make a detailed 
discussion on how these self-consistency conditions are 



a(r, {Ci}) 




FIG. 2: Pictorial representation of the self-consistency of the 
scale-setting method through the universal coupling function 
a(r,{ci}). The point A with a closed path represents the 
operation of Renexivity. The paths BC and CB represent 
the operation of Symmetry, and the paths DF,FE and DE 
represent the operation of Transitivity. 



satisfied or broken by the two frequently adopted scale- 
setting methods: BLM/PMC and PMS. As for FAC, its 
FAC scale is determined by requiring all higher order 
corrections to be zero 3 . FAC satisfies all the above men- 
tioned self-consistenct requirements, whose demonstra- 
tion is similar to that of BLM/PMC and is simpler [40| . 
so we will not repeat it here. 



IV. THE PMC SCALE-SETTING 

The PMC provides the principle underlying BLM 
scale-setting, so if not specially stated, we usually treat 
them on equal footing. 

A. What is PMC ? 

In the original BLM paper [fl| , the physical observable 
is expanded as 



1 + (An f + B) 



(28) 



where is the renormalization scale, the rif term is due 
to the quark vacuum polarization. For clarity, we have 
taken the familiar MS'-scheme. When absorbing all the 



3 This method itself is useful to define an effective coupling con- 
stant for a physical process However it will give wrong re- 
sult when applied to QED processes. The FAC forces all higher 
order corrections to vanish and runs the risk of the better ap- 
proximation being "dragged down" by the poorer one Q. 
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NLO terms involving n/ into the running coupling, we 
obtain fill 



l + C{ 



7T 



where 



33 

fi* = nexp (3A) and C* = — A + B 



(29) 



(30) 



The new scale fj,* and the coefficient Cj" are rif indepen- 
dent. The term 33^4/2 in C* serves to remove that part 
of the constant B which renormalizes the NLO coupling 
constant. 

Through these procedures, it was suggested that the 



pQCD convergence can be greatly improved How- 
ever, after a proper extension of BLM, it can do much 
more than that. 

In deriving Eq. ([2U)) . Brodsky-Lepage-Mackenzie al- 
ready observed that to derive the correct scheme- 
independent LO QED/QCD scale, one should deal with 
the /3o-term rather than the ro^-term. This point has 
lately been emphasized in Refs. [EH [HJ], where an in- 
teresting feature for the NLO BFKL Pomeron intercept 
function u>(Q 2 ,0) has been found; i.e. after using BLM 
scale-setting, the function lj(Q 2 ,0) has a very weak de- 
pendence on the gluon virtuality Q 2 in comparison with 
that derived from the conventional scale-setting under 
the MOM scheme and WS scheme [Hi]. The BLM has 
also been applied with some modifications for determin- 
ing the effective scale in lattice perturbative theory by 
Lepage and Mackenzie (56|, which greatly enhances the 
predictive power of lattice perturbative theory. How- 
ever, BLM in its original form is difficult to be applied 
to higher order calculations because of the emergence of 
higher order n/-terms as n^-term, n^-term, etc.. 

As an extension of BLM scale-setting, a program to 
deal with higher order n^-terms associated with renor- 
malization has been raised in Ref. [36| . which suggests 
that one can expand the effective scale as a perturbative 
series. Later on, an enhanced discussion of this sugges- 
tion up to NNLO level has been presented in Ref. [20j . 
where the perturbative series of the effective scale is ex- 
ponentiated, which is consistent with PMC procedure. 
In that work it is pointed out that the n^-term at the 
NNLO should be first identified with /3g-term and then 
be absorbed into the coupling constant 4 . 

The pioneering work for PMC is done in Ref. 29|, 
which shows that a single global PMC scale, valid at 
LO, can be derived from basic properties of the pertur- 
bative QCD cross-section. Later on, explicit formulae for 



Choose renormalization scheme; e.g. a^f/^t) 
1 

Choose fiinit't arbitrary initial renormalization scale 

1 

Identify {0^} — terms using nf — terms 
through the PMC — BLM correspondence principle 

I 

Shift scale of to ppMC to eliminate } — terms 
Conformal Series 
Result is independent of Ui n i t and scheme at fixed order 



FIG. 3: A "flow chart" which illustrates the PMC procedure, 
where 1Z stands for an arbitrary renormalization scheme. 



scttingPMC scales up to NNLO has been presented in 
Ref. [30] - It has also been pointed out that by intro- 
ducing the PMC-BLM correspondence principle, we can 
improve the previous BLM procedure to deal with the 
process up to all orders, whose estimation is the same 
as PMC. In this sense, we say that PMC and BLM are 
equivalent to each other. Recently, by applying PMC to 
the top-quark pair hadroproduction up to NNLO level 
at the Tevatron and LHC colliders, the most striking 
feature of PMC has been observed, which shows that 
the PMC scales and the resulting finite-order PMC pre- 
dictions are both to high accuracy independent of the 
choice of an initial renormalization scale, consistent with 
RG-invariance [3, H^, [H[ . This implies that the serious 
systematic renormalization scale error introduced by us- 
ing conventional scale-setting can be eliminated by PMC 
through a self-consistency way. 

A "flow chart" which illustrates the PMC procedure 
is presented in Fig. (J3j) , where 1Z stands for an arbitrary 
renormalization scheme. The PMC provides a unambigu- 
ous and systematic way to set the optimized renormal- 
ization scale up to all orders; i.e. we first arrange all the 
coefficients, which usually are given as a series in nj, for 
each perturbative order into {/3^}-terms or non-{/3j^}- 
terms, and absorb different types of {/^j-term into the 
running coupling constant, order-by-order . Different 
types of {/3^}-term are absorbed into different PMC 
scales. Different skeleton graphs can have different PMC 
scales. The PMC scales themselves will be a perturba- 
tive expansion series in a s . After this procedure, all non- 



4 Strictly, together with the n^-term at the same order, it should 
be arranged into a proper linear combination of /3i-term and /3q- 
term, and the /Jg-term will be absorbed into LO PMC scale and 
^i-term will be absorbed into NLO PMC scale pjoll . 



In practice, we can directly deal with ny-terms of the coefficients 
without changing them into {/9?^}-terms, and eliminate the ny- 
terms from the highest power to none also in an order-by-order 
manner. The results are the same due to the PMC-BLM corre- 
spondence [3(ill . 
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conformal {/3^}-terms in the perturbative expansion are 
rcsummed into the running couplings so that the remain- 
ing terms in the perturbative series are identical to that 
of a conformal theory; i.e., the corresponding theory with 
{/3f} = {0}. 

As a simple explanation of PMC, for the coefficient 
Ci(fi) at the NLO level, we have 

Ci(» = C\ {n) +Cn(n)nf, (31) 
= Ci (m)+Cii(m)/3o (32) 

where fi stands for an arbitrary initial renormaliza- 
tion scale, the coefficients Cio(/i) and Cn(/x) are nj- 
independent, C\q = Cio + ^Cn, and Cn = — fCn- The 
LO PMC scale ^pmc is then set by the condition 

Cii(mpmc) = 0. (33) 

This prescription ensures that, as in QED, vacuum po- 
larization contributions due to the light-fcrmion pairs are 
absorbed into the coupling constant. Note that because 
Cxi of Cn, one can practically obtain the PMC scale 
by using the equation Cxx(/-*pmc) = 0, which is usually 
adopted in the literature 6 . However one should keep in 
mind that Eq. (l33l) is exact. 

The PMC - BLM correspondence principle suggested 
in Ref. [3(3] is based on the fact that the purpose of the 
running coupling in any gauge theory is to sum up all 
the terms involving the {/3^}-functions; conversely, one 
can find all the needed {/3^}-terms at any relevant or- 
der by identifying terms arising from the order-by-order 
expansion of the running coupling. This principle pro- 
vides a one-to-one correspondence between the n/-series 
and the /3^-series, and it provide a practical way of iden- 
tifying the terms in the n^-series into the required 0^- 
series. The {/3^}-series derived from Eq. (fT3"|) provides 
the foundation for the PMC - BLM correspondence prin- 
ciple, since it shows which {/^j-terms should be kept at 
a specific perturbative order. This procedure provides a 
convenient and consistent way of treating the {/3^}-terms 
in the perturbative series. Its advantages will be shown in 
the next subsection. Such a choice of {/^j-series is not 
completely identical to the suggestion of Refs. [36h38u57| . 
In Refs. 36-38], as an extension of BLM scale-setting to 
all orders, the large /3o-approximation is adopted with 
some modifications to simplify the calculation (called as 
the seBLM and the xBLM approach |37j) 7 . 



6 This should be used with care, since if Cio is a constant free of 
scale, then such practical way will give wrong NLO coefficient 
other than the correct one Cio. 

7 Theoretical differences for different treatments will be discussed 
in more detail and will be presented elsewhere. 



B. The properties of PMC 

It is straightforward to verify that PMC satisfies all 
the self-consistency requirements outlined above. 

1. The Existence and Uniqueness of the renor- 
malization scale /x are guaranteed, since the scale- 
setting conditions for PMC are often linear equa- 
tions in In /it 2 . 

As a simpler explanation, if the NLO coefficient 
Ci (fi) in Eq. (fTJ) has the form 

Ci(fi) = (a + b n f ) + (c + d n f )\n fi 2 , (34) 

with a, b, c and d are constants free of n/ . The LO 
PMC scale can be set as 

hx^ c = ~ + 0{a s ), (35) 

where the higher order a s -terms will be determined 
by n/-terms at the NLO-level or even higher levels. 

2. Refiexivity is satisfied. The PMC requires all 
ln(/x 2 //x /2 )-terms in Eq. (Tl4]) vanish, thus we obtain 

fi' = fj, . 

3. Symmetry is trivial, because after PMC scale- 
setting, we always have 

ru(in,i4) = -hid**, i4) ■ 

That is, the two NLO coefficients only differ by 
a sign. Thus, requiring one of them to be 
independent is equivalent to requiring the other one 
also to be {/J^j-independent. This argument en- 
sures the symmetric relation, A12A21 = 1, is satis- 
fied after PMC scale-setting. 

4. Transitivity is also satisfied by PMC. After PMC 
scale-setting, the two coefficients f 12(111, //J) and 
^23(^21 M3) m t ne following two series 

= a a2 {ti)+h2{vi,ti)a 2 M + 0{al 2 ) (36) 
and 

0^2(^2) 

= a s3 (n* 3 ) + f 23 (^, V* 3 )a 2 s3 (fi* 3 ) + 0(e&) . (37) 

should be independent of {Pi}. After substituting 
Ea.([37]) into Eq.®, we obtain 

= Ot s3 (f4) + [f 12^1, f4) + ^23(M2>M3)] Q 4(A t 3) 

+<D(a 3 s3 ) . (38) 

We see that the new NLO coefficient [fi 2 (^i, + 
^23(^*2% M3)] wn l a lso be {/^j- independent, since 
it is the sum of two {/^j-independent quanti- 
ties. These arguments ensure the transitive rela- 
tion, A31 = A32A21, be satisfied after PMC scale- 
setting. 
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As a combination of all the above mentioned PMC 
features, the advantages of PMC are clear 8 : 

• It keeps the information of the higher order cor- 
rections but in a more convergent perturbative se- 
ries. After PMC scale-setting, the divergent "renor- 
malon" series with n!-growth disappear in the per- 
turbative series, so that a more convergent pertur- 
bative series is obtained. Such better convergence 
has already been found in the original BLM pa- 
per [lH and the following BLM-literature even at 
the NLO level. 

• After PMC scale-setting, the renormalization scale 
dependence is transformed to the initial renor- 
malization scale dependence, and it is found that 
such initial renormalization scale dependence can 
be highly suppressed or even eliminated : 

— The resulting expressions are conformally in- 
variant and thus do not depend on the choice 
of renormalization scheme. 

— One can obtain proper scale-displacements 
among the PMC scales which are derived un- 
der different schemes or conventions. 

— One can obtain the CSR between any two 
physical observables such as the Generalized 
Crewther Relation connecting the Bjorkcn 
sum rule to the e + e~ annihilation cross sec- 
tion. Many leading order CSRs have been de- 
rived in Ref. [13] ■ The CSRs have no scale am- 
biguity and are independent of the choice of 
renormalization scheme. The relative scales 
in the CSR ensure that two observables pass 
through each quark threshold in synchrony. 
The coefficients in the CSR can be identified 
with those obtained in conformally invariant 
gauge theory (H, HHIH. 

— There can be some residual scheme depen- 
dence for a fixed-order calculation due to 
unknown higher-order terms. However this 
scheme-dependence can be highly suppressed 
in a similar way as that of the residual ini- 
tial scale dependence; such effects can be es- 
timated by using the RG-based scheme equa- 
tions [13, 53. 

• The PMC provides a fundamental and system- 
atic way to set the optimized renormalization scale 



In the PMC, the same procedure is valid for both space-like and 
time-like arguments; in particular this leads to well-behaved per- 
turbative expansion, since all the large {/3^}-dependent terms on 
the time-like side involving 7r 2 -terms are fully absorbed into the 
coupling. The PMC does not change the space-like or time-like 
nature of the initial renormalization scale Qo , since in general, all 
the PMC scales are equal to Qo times an exponential factor |3(J1 . 



for the fixed-order calculation. In principle, PMC 
needs an initial renormalization scale to initialize 
it. However, it is found that the estimates after 
PMC scale-setting are independent of any choice of 
the initial renormalization scale - even the PMC 
scales themselves are independent of any choice 
of initial scale and are 'physical' at any fixed or- 
der. This is because that the PMC scale itself is a 
perturbative series and the unknown higher-order 
{/3^}-terms are to be absorbed into the higher- 
order term of PMC scale and will be strongly 
power suppressed. One example of this behav- 
ior is shown in Refs. [H], [H, H3, where the top- 
quark pair total cross-section and the top-quark 
pair forward-backward asymmetry are almost free 
from the choice of initial renormalization scale even 
at the NNLO-level. 

• Moreover, it is found that the PMC scale-setting 
can also be adopted for QED case. The variable 
Nc can be taken as an analytic variable. In the 
Abelian limit Nc — >• at fixed a = Cf&s with 
C F = (7V C 2 - l)/2N c Hi, the PMC method agrees 
with the standard Gell Mann-Low procedure for 
setting the renormalization scale in QED, a con- 
sistency requirement of analyticity of Yang Mills 
gauge theories. 

• After PMC scale-setting, the number of active fla- 
vors rif is correctly determined [48| . Using the 
PMC ensures that the expansion is unchanged as 
one passes each quark threshold, since all vacuum- 
polarization effects due to each new quark are auto- 
matically absorbed into the effective coupling con- 
stant. 

• The argument of the running coupling has time- 
like or space-like values appropriate to the physics 
of the PMC scale; for example the scale of the QED 
coupling which some all vacuum polarization cor- 
rections in the lowest order e + e~ — > ampli- 
tude is a(t) in the Gell Mann-Low scheme. As in 
QED, the running QCD coupling is complex in the 
time-like domain, reflecting the contribution of di- 
agrams with physical unitarity cuts. 



V. THE PMS SCALE-SETTING 
A. What is PMS ? 

The PMS states that 0413] if an estimate depends on 
some "unphysical" parameters 9 , then their values should 
be chosen so as to minimize the sensitivity of the estimate 



Here the "unphysical" parameter means which is known not to 
affect the true result. 
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to small variations in these parameters; i.e. this method 
chooses /ipMS a t the stationary point of pjv : 



Op 



N 



dpN i 







d\a(p?) 



. 



(39) 



(40) 



Here Eq. (|4"0|) can be solved with the help of the usual 
renormalization group equation (j4|. 



B. The properties of PMS 

Unlike the case of PMC, in general, there are no known 
theorems that guarantee the Existence or the Unique- 
ness of the PMS solution. Although for practical cases 
PMS does provide solutions, and when there are more 
than one solution usually only one of them lies in the 
physically reasonable region these observations 

alone do not guarantee that PMS will be trouble-free for 
new processes. 

To discuss PMS properties in a renormalization 
scheme-independent way, following the suggestion of 
Ref. [4(j, we adopt the 't Hooft scheme 62 1 to define 
the running behavior of the effective coupling constant. 
Under the 't Hooft scheme, all the scheme parameters 
{cj} are set to zero, and Eq.© simplifies to 



da 



-a 2 (l + a), 



whose solution can be written as 



1 



t = — h In , 
a \ 1 



(41) 



(42) 



In the above solution, for convenience, we have redefined 
where ^A^d ' 1S the asymptotic scale 



t as hi 



A tH2 

QCD 



under the 't Hooft scheme. The 't Hooft coupling con- 
stant has a formal singularity, a(r, {c,}) = a(0, {0}) = 
co, which provides a precise definition for the asymp- 
totic scale Aq£ D (62[; i.e. it is defined to be the pole of 
the coupling function. 

Given two effective coupling constants a\ and <Z2 under 
the 't Hooft scheme, they are related by the perturbative 
series 



ai(ri) = a 2 (r 2 ) + (t 2 - n)ai( T 2) + 



(43) 



PMS proposes the choice of \ii (or equivalently, r 2 ) at 
the stationary point, i.e.: 

P- = = [03(75) + (t-2 - 71)^(75)] . (44) 
«t 2 ar 2 

Then, we obtain the condition: 



In order to obtain 72 in terms of ri, one must solve the 
last equation in conjunction with 



1 

«2 



log 



a 2 



l + a 2 



r 2 . 



(46) 




FIG. 4: The dependence of the PMS scale parameter T2 as a 
function the external scale parameter n. 

In Fig[4]we present the graphical solution of the PMS 
scale-parameter 72 as a function of the external scale- 
parameter t\. One may observe two points: 

• 72 > 7i — 5. Since r 2 ^ t\ in any cases, so PMS 
explicitly violates the Reflexivity. For a fixed- 
order estimation, when one uses an effective cou- 
pling constant to predict itself, the application of 
PMS would lead to an inaccurate result. 

• In the large momentum region (n S> 1), we obtain 
02(7-2) -¥ 0, and 



72 



1 



(47) 



Under the same renormalization scheme R, we have 
the same asymptotic parameter ^-qcd R f° r both 

cii and a,2- Here Ag^^ fl is the 't Hooft scale as- 
sociated with the i?-scheme, where the word "as- 
sociated" means we are choosing the particular 't 
Hooft scheme that shares the same 't Hooft scale 
with the i?-scheme. Then the relation (j47j) in terms 
of Hi and \ii becomes 



fi 2 ~ fii exp 



(48) 



More generally, it is found that after PMS scale- 
setting, the scale displacement between any two scales 
Hi and pj in the large momentum region is 
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This would mean that 

A 12 A 21 ~ exp (~^) ^ 1, (50) 

A13A32A21 ~ exp ^ L ( 51 ) 

This shows that the PMS does not satisfy the Symme- 
try and Transitivity requirements. Let us point out 
that adding the scheme-parameter optimization in PMS 
does not change any of the above conclusions. It only 
makes the solution much more complicated [631 ] - The 
inability of PMS to meet these self-consistency require- 
ments resides in that the derivative operations in gen- 
eral do not commute with the operations of Reflexivity, 
Symmetry and Transitivity. 

As argued in Sec. Ill, any truncated perturbative series 
will explicitly break RG-invariance i.e. Eq.© can 
only be approximately satisfied for any fixed-order esti- 
mation. The precision depends on to which perturbative 
order we have calculated, the convergence of the pertur- 
bative series, and how we set the renormalization scale. 
As shown by Eq. pp|) . the PMS requires the truncated 
series, i.e. the approximant of a physical observable, to 
satisfy the RG-invariance near /i = /ipMS- This provides 
the underlying reason for why PMS docs not satisfy the 
Reflexivity, Symmetry and Transitivity properties. 
This shows the necessity of further careful studies of the- 
oretical principles lying below PMS. 

The PMC and PMS scale-setting methods each gives 
specific predictions for physical observables at finite or- 
der; however, their predictions are very different: 

• The PMC sums all {/3^}-terms in an arbitrary 
renormalization scheme 1Z in the fixed-order pre- 
diction into the running coupling, leaving the con- 
formal series. It satisfies all of the RG-properties 
Reflexivity, Symmetry, and Transitivity. The 
PMC prediction is thus scheme- independent, and 
it automatically assigns the correct displacement 
of the intrinsic scales between schemes. The vari- 
ation of the prediction away from the PMC scale 
exposes the non-zero {/^j-dependent terms. The 
PMC prediction does have small residual depen- 
dence on the initial choice of scale due to the trun- 
cated unknown higher order {/3^}-terms, which 
will be highly suppressed by proper choice of PMC 
scales. 

• The PMS chooses the renormalization scale such 
that the first derivative of the fixed-order calcula- 
tion with respect to the scale vanishes. This crite- 
rion of minimal sensitivity gives predictions which 
are not the same as the conformal prediction, and 
the PMS prediction depends on the choice of renor- 
malization scheme 10 , and it disagrees with QED 



scale-setting in the Abelian limit. For example, in 
the case of e + e~ — > gqq, the PMS scale decreases 
with increasing gluon jet mass and increasing fla- 
vor number, opposite to the correct physical be- 
havior [si]. The PMS does not satisfy the RG- 
properties of Symmetry, Reflexivity, and Tran- 
sitivity, so that relations between observables de- 
pend on the choice of the intermediate renormal- 
ization scheme. 



VI. SUMMARY 

The conventional scale-setting procedure assigns an ar- 
bitrary range and an arbitrary systematic error to fixed- 
order pQCD predictions. As we have discussed in this ar- 
ticle, this ad hoc assignment of the range and associated 
systematic error is unnecessary and can be eliminated by 
a proper scale-setting method. 

Renormalization group invariance ([9]) states that a 
physical quantity should be independent of the renormal- 
ization scale and renormalization scheme. With the help 
of the RG based equations which incorporate the scheme 
parameters, we have presented a general demonstration 
for the RG-invariance by setting the perturbative series 
up to infinite orders. 

We have discussed the necessary self-consistency con- 
ditions for a scale-setting method, such as the Existence 
and Uniqueness of the renormalization scale, Reflex- 
ivity, Symmetry, and Transitivity. There properties 
are natural deductions of RG-invariance. We have shown 
that PMC satisfies these requirements, whereas the PMS 
does not. We have also pictorially argued that the formal- 
ism based on the RG based equations satisfies all these 
requirements for scale and scheme variation. 

The principle of minimum sensitivity (PMS) requires 
that the slope of the approximant of an observable to 
vanish at the renormalization point. With the help of the 
RG based equations, it has been argued that PMS can 
provide renormalization-scheme dependent estimates 0- 
[l(j. We have shown that the PMS violates the Sym- 
metry and Transitivity properties of the renormaliza- 
tion group, and it does not reproduce the Gell-Mann- 
Low scale for QED observables. Eq. (|5"Tj) shows that the 
relation between any two physical observables after PMS 
scale-setting depends on which renormalization scheme 
chosen for the calculation, which explicitly breaks the 
"group properties" of the RG equations. In addition, the 
application of PMS to jet production gives unphysical re- 
sults [64| , since it sums physics into the running coupling 
not associated with renormalization. This implies the ne- 
cessity of further careful studies of theoretical principles 



As shown in Ref. 1631 . by using the PMS together with the 



scheme-equations J7J and the scheme-independent equation 1101 . 
such renormalization scheme dependence can be reduced to a cer- 
tain degree through an order-by-order procedure. 
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lying below PMS. 

The principle of maximum conformality (PMC) pro- 
vides a fundamental and systematic way to set the op- 
timized renormalization scale at fixed order in pQCD. 
The PMC has a solid theoretical background |29|, |30| . 
it provides the underlying principle for BLM, and many 
PMC features have already been noted in the BLM liter- 
ature. Most important, it is found after standard PMC 
scale-setting, the theoretical prediction is essentially in- 
dependent of the choice of initial renormalization scale 
and the theorist's choice of renormalization scheme, con- 
sistent with the RG invariance. 

The most important goal for a scale-setting method is 
to eliminate the renormalization scheme and initial scale 
dependences - more fundamental requirements than im- 
proving convergence of the pQCD series. In the liter- 
ature, however, some extensions of BLM scale-setting 
have concentrated on how to improve the pQCD conver- 
gence, such as the large /^-expansion [65l |. the sequential 
BLM (seBLM) and xBLM [36M38I] . etc.. In fact, once 
one sets the scales properly, as PMC does, much better 
pQCD convergence than the conventional scale-setting 
method is automatic, since the divergent "renormalon" 
series with n!-growth has been absorbed into the effective 
scales and disappears in the perturbative series. An ex- 
ample of this improved convergence can be found in our 
analysis for the top-quark pair production at the NNLO 
level [13, [H, HH . 

Two more subtle points for PMC scale-setting: 

• In some specific kinematical regions, such as for 
the heavy quark pair production in the threshold 
region, Coulomb-type corrections will lead to siz- 
able contributions which are enhanced by factors 
of ir/v and the PMC scale can be relatively soft 
for v — ¥ (v, the heavy quark velocity). Thus the 
terms which are proportional to (tt/v) or (ir/v) 2 
should be treated separately in that different PMC 
scales are adopted in the estimation (2|| [66| . 

• The factorization scale fi f which enters into the pre- 



dictions for QCD inclusive reactions is introduced 
to match nonperturbative and perturbative aspects 
of the parton distributions in hadrons. The fac- 
torization scale occurs even for a conformal theory 
with {fij 1 } — where a s is constant. The factor- 
ization scale should be chosen to match the non- 
perturbative bound state dynamics with perturba- 
tive DGLAP evolution. This can be done explicitly 
for electron-atom or atom-atom inelastic scatter- 
ing processes in QED using the known bound state 
dynamics of atoms. This could also be done in 
hadron physics using nonperturbative models such 
as AdS/QCD and light- front h olog raphy; recent re- 
views can be found in Refs. [o7l.l68l] There is clearly 
no reason to equate the factorization scale to the 
renormalization scale [69[ . We expect that the fac- 
torization scale ambiguity will also be reduced by 
applying the PMC scale-setting to the kernels of 
DGLAP evolution equations. 

In summary, the systematic application of the PMC 
can eliminate a major ambiguity of pQCD predictions 
from scale and scheme ambiguities, thus greatly improv- 
ing the precision of tests of the Standard Model and the 
sensitivity to new physics at colliders and other experi- 
ments. 
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